Low-Lying Excitations of Quantum Hall Droplets 
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We present the adequate circular representation of charge 
density excitations in Quantum Hall droplets at filling factor 
1. A novel set of operators describes magnetoexcitons with 
well defined both total and center-of-mass angular momen- 
tum. The accuracy of this description is shown by the high 
overlap of the magnetoexcitons with the exact eigenstates of 
a quantum dot in a magnetic field obtained from numerical 
calculations up to 20 electrons. From such magnetoexcitons 
we get a good understanding of total energies and spectral 
functions. 

PACS number: 73.40.Hm 

The electronic structure of quantum Hall droplets 
(QHD) has been studied extensively |[l|-p^, both due 
to its relevance in the analysis of quasi-two dimensional 
(2D) artificial atoms and in order to get insight 

on the behaviour of 2D systems in the quantum Hall 
regime. Special attention has been given to 2D parabolic 
quantum dots (PQD) since the non-interacting system 
is analytically solvable, its high symmetry symplifies the 
calculations and some of its properties do not depend on 
the strengh of the confinement potential. 

In this paper we introduce the set of operators describ- 
ing low- lying excitations in PQD for filling factor v ~ I. 
There are two key conditions to obtain these operators: 
Firstly, the fact that, close to v = 1, the lowest excita- 
tion within each total angular momentum subspace are 
excitations of integer charge and, secondly, to impose to 
the eigenstates these operators generate a well defined 
third component of both total (M) and center-of-mass 
{Mcm) angular momentum. At the end of the paper we 
discuss several applications of such operators in the anal- 
ysis and calculation of relevant magnitudes for transport 
and spectroscopy in QHD. 

We consider interacting electrons in a PQD affected by 
a perpendicular magnetic field B high enough for the N- 
electron wave functions to comprise only polarised single 
particle states in the lowest Landau level (LLL). Due to 
the circular symmetry, it is convenient to work in the 
symmetric gauge, so that M is a good quantum number. 
The Hamiltonian can be written as: 
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where the sum runs over the single-particle angular mo- 
mentum labelling the states of a PQD defined in the 
xy plane, and characterized by a frequency wq- 7 = 



uj ^/Q., = { eB)/{m*c), Vl = {uj^ +4^2)1/2^ and Ig = 
\/ {h)/{m*n) is the effective magnetic length, c' and c 
are the electron creation and annihilation operators. An 
additional advantage of PQD is that the electron-electron 
interaction matrix elements Vmim2mami between states 
{to} in 2D can be computed analytically |p^,[7|. In this 
paper we take the interaction as coulombic, unless ex- 
plicitely stated otherwise. Let us review some of the 
known properties of this Hamiltonian. Since the inter- 
action only depends on the relative motion of electrons, 
the QHD spectrum separates in subspaces with differ- 
ent value of M [^,^,|l^. Due to the reduction of the 
Hilbert space to the LLL, the dimension of each sub- 
space labelled by {N, M) is finite. Remarkably, wave 
functions do not depend on wq or B (apart from a rescal- 
ing of the length scale), neither does their energy ordering 
within a {N, M) subspace. This property can be used to 
get information on an infinite system from the results 
in a PQD. Only the relative positioning of energies for 
different [N, M) subspaces depend on the values of the 
parameters defining the QD due to the first term of the 
Hamiltonian. On increasing B (or, alternatively, decreas- 
ing Wo) the A^-electron GS resides in {N, M) subspaces 
with increasingly larger M . 



We denote the GS in subspace {N,M) by | $ 
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An important fact, which to the best of our knowledge, 
has never been stated before, is that {Ne)ai = {^m^ \ 
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is extremely close to an integer. 



As M increases, {Ne)m increases in steps of heigh 1, 
and size typically N, up to a number of steps of the or- 
der of N/2 beyond which this property fails jl|]. We 
have checked that this property holds not only for the 
Coulomb interaction, but also for the logarithmic and 
contact interactions. Therefore, in order to study the 
low-lying excitations close to Mq, we can consider a fam- 
ily of variational wave functions, each one of them being a 
linear combination of Slater determinants with the same 
value of {Ne)m- This additional "symmetry" allows the 
computation of properties in this system for very large 
values of N and M. 

In this letter we are going to concentrate on the study 
of the first sector, i.e. when {Ne)m ~ 1- Numerically, we 
find that there is always an eigenstate with the previous 
property in each subspace {N, M) with Mq < M < Mq + 
N. Excitations with {Ne)m = 1 can be expressed as a 
representation of the charge density operator projected 
onto the LLL. They resemble the single-mode excitations 
previously studied in superfluid *He and in the fractional 
Quantum Hall regime Mm . A set of such operators which 



has been put forward [[4y5|,|7|,|ll|-|13| to describe boson-like 
edge- wave excitations in PQD is 



N 
k=l 



h\ = -(xfc -I- iyk)/'2.l 



B — —Zk/2 is, apart from a normal- 
ization constant, the ladder operator increasing the single 
particle angular momentum in the LLL. It has been con- 
jetured 0, and proven for the parabolic electron-electron 
interaction , that these operators represent a bosoni- 
sation of the low- lying excitations for AM « N . How- 
ever the overlap of the states they generate with the true 
ground state decreases rapidly with AM being negligible 
for AM K, N (see below). The reason for this is that 



•^AM I *^Mo^)' ^'^^ AAf > 1, is not an eigenstate of the 
third component of the center- of-mass angular momen- 
tum. Note that S'|, when operated on an exact eigen- 
state, gives an eigenstate with Mcm increased by 1 [||. 
States with well defined Mcm should be anihilated by 
(^S^'^McM _ ji; jg straighforward to check that the condition 

(^Si)McM sf^^^ I $2,^) = is not fulfilled for AM > 1. 

The previous discussion gives the key idea for finding 
boson-like excitations (with well defined M and Mcm) 
by means of an operator J^j^j which ejects one electron 
from the compact state. Since these excitations must 
have the lowest possible energy, the center-of-mass tends 
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to be at rest, implying the condition SiJ^^^j \ -±-^^^ 
From this condition we get, 
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for 1 < AM < N . These operators can be written as 

Am 2'^^'^/2/dz(z*)-'^^Vo(2)- They are the repre- 
sentation of the projection on the LLL of the density 
operator (/9o(-z)) on the subset ((z*)"^*'^) of the cylin- 
drical harmonics |po|| . From the first quantisation ex- 
pression for J^^,/, we find ||2l| that it is possible to ex- 
pand J^jvf in a series which i) for AM = >> 1 
gives the Laughlin [0 wavefunction for a quasi-electron: 

I ^^N^) ~ ri-^^j I ^Mo^)' s-'^d ii) has S\j^i as asymtotic 
limit for AM << N. In between these two limits the 
density charge excitation created by the operator J^^j 
consists of a "dressed" electron- hole pair, i.e., a magne- 
toexciton, where the distance between the hole and the 
electron is proportional to the quantum number AAI. 
The electron is always located at the edge of the droplet 
while the hole moves toward the center of the droplet 
with AM. For AAI — N, the hole reaches the centre of 
the droplet and this marks the end of the first sector. 

In order to check the quality of J^j^.f for de- 
scribing QHD excitations, table I shows the overlaps 
I^ImI<^) and |5L,|<)) for N = 10,20 

and AM < N. Note that for AM < 3, JImI^Mo^) is 



an exact eigenstate. Using Trugman and Kivclson's Q| 
procedure to generate eigenstates it is possible to show 
that this is so for any isotropic pairwise interaction. The 
conclusion from table I is clear: the branch of lowest 
excitations are magnetoexcitons with the center-of-mass 
at rest. S'l a/ only describe approximately the excita- 



tions in the AM < N limit |5|j3,|ll|-|13f and, as expected, 
become poor descriptions of the excitations for droplets 
with small number of electrons. The operators J^^/ rep- 
resent a bosonization of low-lying excitations of the PQD 
close to = 1 . As it happens in other cases of bosoniza- 
tion they do not satisfy bosonic commutation rela- 
tions strictly; however their expectation values on | '1'^/'^) 
do. ' 

From the analytical expression for the wavefunc- 



tions it is possible to compute their energies E 
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large values of N. Using Wick's theorem, the energy 
difference between the GS in the {N,M) and {N,Mo) 
subspaces can be written as: 
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for 



nn{l 7)AM/2 + 4^) (AM) + 6n(AA^) + C;.(AM). 

The first term is the single-particle energy difference be- 
tween the two states. The remaining three terms are the 
contributions coming from the interaction; they involved 
combinations of interaction matrix elements, and can be 
identified as: e^^^ is the interaction between the com- 
pact V = 1 droplet and the ejected electron onto the edge 
and is practically independent on AM. e^^)j is the en- 
ergy corresponding to the interaction between the v = 1 
droplet and the hole created in it, and decreases almost 
linearly with AM. tf^f^ represents the binding energy of 
the magnetoexciton, and presents a minimum as a func- 
tion of M. Figure |l] shows the dependence with M of the 
these last three terms for N — 30. We find that the total 
energy difference has the magnetoroton minimum at the 
same position as Cg^j, in the region where the description 
in terms of valid. The magnetoroton minimum 

marks the AM at which the first instability in the v = 1 
droplet occurs. We find that, for a contact interaction, 
the position of this intra-Landau magnetoroton at v = 1 
shifts as 2\/]V in the whole range of N. 

We consider now another important application of the 
operators Ja j\,t- Several transport and spectroscopic 
properties fi p3| , p^ are controlled, close to = 1, by 
the spectral function 
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where rj labels all A^-particle states. In spite of the fact 
that magnetoexcitons •/amI^^Mq^) ^I'e not the only exci- 
tations in the system, the main contributions to A[N,uj) 

come from states of the form (<S'J)-' J^m I ^'mq)^ thai is, 
from magnetoexcitons and all the states generated from 



them by increasing Mcm- Figure ^ shows the compari- 
son between the exact spectral function A{N,lo) (lower 
panel) and the one calculated including only the preceed- 
ing states (upper panel), for = 8 and at a value of 
for which the GS of the 7-electron PQD is the compact 
state. The identification of the relevant states for the cal- 
culation of A{N, Lo) allows i) the analytical calculation of 
this magnitud for large N and ii) the understanding of 
the different bunches of peaks spaced in (Jl — lOc)/2 and 
decreasing in intensity that appear in A{N,(jj). After the 
first high peak coming from the transition between com- 
pact states in both A^— 1 and TV-particle systems, the first 
bunch is due to the transition from the compact state of 

— 1 particles to magnetoexcitons of N particles. Each 
subsequent bunch of peaks in A{N,u)) is a replica of the 
preceeding one with Mcm increased in 1. 

Let us now pay attention to the spectrum for large val- 
ues of AAf. The lowest-lying states of sectors beyond the 
first, and even at the end of the first sector (AM « N) for 
large values of N , cannot be described in terms of single 
magnetoexcitons. This is observed in table I where the 
overlap at AM = TV = 20 is small Since an integer 
number of electrons is ejected from the compact state, 
one is tempted to look for almost independent bosons in 
which low-lying states are obtained by succesive applica- 
tion of J^jvf- Comparing the energies so obtained with 
the exact ones, the result is satisfactory. However this 
is not generally the case when comparing wavefunctions. 

The state Jaj\/i -^aa/s I*^Mo^) ^i*'^ ^-^i + ^^^2 > N has 
a small overlap with the state |*I'm1+aMi+am2) obtained 
by exact diagonalization. Only when AMi = AM2 = N , 
the overlap is close to 1 and the state can be considered 
as constituted by two almost non-interacting magnetoex- 
citons. 

Finally, let us discuss one implication of the increas- 
ing correlation effects with increasing B. GS changes 
from the compact state to a magnetoexciton and fur- 
ther to states with several holes in the compact and 
the same number of electrons out of it. Once spectral 
functions are known, single electron capacitance exper- 
iments in QD can be analized. For low tempera- 
ture, tunneling rates depend only on the spectral weight 
A(iV_) = A{N,^i{N)) 0. This magnitude is shown in 
fig. 1^ as a function of B. For a field in which the GS 
of both TV — 1 and N electron systems are the compact 
ones, the spectral weight is obviously 1. For increasing 
field those two GS melt by ejecting 1 electron. As holes 
produced in the inner part of the dot are not completely 
equivalent, A(iV) decreases. In some cases the value of 
the field for which melting appears is significantly differ- 
ent between iV — 1 and iV particles so that the spectral 
weight is drastically quenched. The reduction of A(Af) 
with B explains the experimental quenching of signal 
amplitude in single electron capacitance experiments for 
high magnetic fields fl^ . 

In summary, we have found the adequate circular 
representation of charge density excitations in QHD in 



u = \. Magnetoexcitons J ]\m\^ Mo) turn out to be an 
excellent description of low- lying excitations. From such 
wavefunctions it is easy to get energies and spectral func- 
tions, i. e. the important information on the electronic 
structure of QHD. 
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FIG. 1. Different contributions to the magnetoexciton en- 
ergy as a function of AM for 30 particles. 



FIG. 2. (a) Spectral function for 8 electrons includ- 
ing only magnetoexcitons with rotating center-of-mass 

((5'J)-' J^jyj|<I>^^')). (b) Total spectral function including the 
whole spectrum obtained from the exact diagonalization. 

FIG. 3. Many-body contribution to the low-temperature 

tunneling rates of an electron entering in a QD as given by the 
spectral weight A (A/") as a function of B for different number 
of electrons. 



